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In this work, we have systematically studied structural, electronic and magnetic properties of 
atomic scale defects in 2D transition metal dichalcogenides MX 2 , (M = Mo and W; X = S, Se and 
Te) by density functional theory. Various types of defects, e.g., X vacancy, X interstitial, M vacancy, 

M interstitial, MX and XX double vacancies have been considered. It has been found that the X 
interstitial has the lowest formation energy (~ 1 eV) for all the systems in the X-rich condition 
whereas for M-rich condition, X vacancy has the lowest formation energy except for MTe 2 systems. 

Both these defects have very high equilibrium defect concentrations at growth temperatures (lOOOK- 
1200K) reported in literature. A pair of defects, e.g., two X vacancies or one M and one X vacancies 
tend to occupy the nearest possible distance. No trace of magnetism has been found for any one of 
the defects considered. Apart from X interstitial, all other defects have defect states appearing in 
the band gap, which can greatly affect the electronic and optical properties of the pristine systems. 

Our calculated optical properties show that the defect states cause optical transitions at ~ 1.0 eV, 
which can be beneficial for light emitting devices. The results of our systematic study are expected 
to guide the experimental nanoengineering of defects to achieve suitable properties related to band 
gap modifications and characterization of defect fingerprints via optical absorption measurements. 


I. INTRODUCTION 

Ever since the first successful fabrication of graphene, 
a single layered honeycomb structure of carbon atoms 
by Novoselov et al in 2004 ^, there has been a contin¬ 
uous increase in research interests on many other 2 D 
monolayers. Apart from graphen^^ other extensively stud¬ 
ied sin gle lay ered systems ye sil icen germanen^^®^^, 
stanene^l^E®, phosphoren d^^ * ^^ * and a number of single 
layered transition metal dichalcogenides (MX2, M = 
Transition metal and X = Chalcogen)^^®^. These single 
layered systems, having two dimensional ( 2 D) structures, 
generally exhibit properties, which are remarkably differ¬ 
ent from their respective three dimensional bulk phases 
and thus these systems have become important subjects 
for experimental and theoretical studies. From the appli¬ 
cation point of view, the 2 D semiconducting systems have 
the minimum possible thickness and wide band gaps, 
which make them important for photovoltaics, sensing, 
biomedicine, solar cells and catalysis. 

It is not uncommon to find atomic scale defects in 
2 D materials, e.g., in graphene, mono and di-vacancies, 
Stone-Wales defects etc. have been studied quite exten¬ 
sively. Vacancies are formed during the fabrication of 
monolayer M0S2 using sonochemical deposition method. 
It is known that these vacancies have considerable effects 
on electronic, magnetic and optical properties. For ex¬ 
ample M0S2 has been found to acquire S vacancies during 
its fabrication. These S vacancies are found to be detri¬ 
mental for n-type conductivi ty of M0S2 as they create 
deep trap states for electron d^^ l ^^ l Moreover, the con¬ 
centration of defects depends on the mode of the fab¬ 
rication process. For example, the carrier mobility in 
M0S2 monolayer fabricated by chemical vapor deposition 
method is 0.02 cm^/V-^^ whereas in mechanically exfo¬ 


liated monolayer, it can go up to 10 cm^/V-^. Owing 
to the importance of role of defects, there are few recent 
theoretical studies of native defects in Ton- 

gay et al have investigated the effect of anion vacancy 
on photoluminescence of M0S2, MoSe2 and WSe2 using 
experiment and density functional theory calculations.^ 
Recently Liu and coworkers have studied the influence of 
anion and metal vacancies on electronic and optical prop¬ 
erties of monolayer MoS2p^ Various defects in monolayer 
M0S2 have been investigated previously by mean s of fi rst 
principles DFT calculations by different group^^^^^^H^. 
Rotational defects have also been studied theoretically 
and experimentally recently on MoSe2 and WSe2 by Lin 
et al. ^ 


However, to the best of our knowledge, systematic 
studies on the role of various point defects and double de¬ 
fects on the electronic and optical properties of transition 
metal dichalcogenides are not available in the literature. 


Therefore, in this paper we have systematically inves¬ 
tigated the role of defects in modifying the electronic 
and optical properties of monolayers of transition metal 
dichalcogenides, MX2 (M = Mo, W and X = S, Se, 
Te). We present the most probable defects under various 
growth conditions, the electronic structures of defected 
materials and the signature of these defects in optical 
properties. The paper is organized in the following way. 
First, we present the computational details followed by 
the results of optimized structures in presence of defects. 
Then we report formation energies and equilibrium defect 
concentrations. Finally, we present electronic structures 
of defected monolayers and their optical signatures. 
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II. COMPUTATIONAL DETAILS 


B. Defect formation energy 


A. Basic parameters 


The defect formation energy is defined as follows. 


All the calculations have been performed with a mono- 
layer MX2 supercell, where M stands for Mo, W and X 
stands for S, Se, Te. We have considered both M and 
X defects (in terms of vacancy or interstitial) in all of 
the above systems. The supercell is generated by repeat¬ 
ing the primitive cell by five times in a and b directions. 
A vacuum of 20 A is included to avoid the interaction 
between periodic images in the out-of-plane direction. 
The calculations have been performed using a plane-wave 
based density functional code VASP.^^ The generalized 
gra dient approximation of Per dew, Burke and Ernzer- 
Pof has been used for the exchange-correlation po¬ 
tential. The structures have been optimized using the 
conjugate gradient method with the forces calculated 
using the Heilman-Feynman theorem. The energy and 
the Heilman-Feynman force thresholds have been kept at 
10“^ eV and 10“^ eV/A respectively. For geometry op¬ 
timizations, a 3 X 3 X 1 Monkhorst-Pack k-giid is used. 
Total energies and electronic structures are calculated 
with the optimized structures on a 5 x 5 x 1 Monkhorst- 
Pack k-grid. It is worth mentioning that some defect 
states may result in finite magnetic moments. For that 
purpose, we have considered spin-polarization in all our 
calculations. However, all the systems relaxed to non 
magnetic ground states. Spin-orbit coupling was not in¬ 
cluded. 


Table I. Comparison of in-plane lattice constant a (in A) and 
band gaps (in eV) for different monolayers of MX 2 as 
calculated with PBF functionals together with exprimental 
values. Experimental values of lattice constants have been 
taken from Refs. E21 and 1331 Experimental values of band 
gaps have been taken from Refs. l34H40l 



M 0 S 2 MoSe 2 MoTe 2 

WS 2 WSe 2 WTe 2 

a (A) PBE 

3.18 3.32 3.56 

3.18 3.32 3.56 

Expt. 

3.16 3.30 3.52 

3.15 3.28 

Eg{eV) PBE 

1.73 1.47 1.08 

1.89 1.58 1.09 

Expt. 

1.95 1.58 1.1 

1.99 1.65 


Ef — F/defect \_E] 


pristine 


, ni/ii 


( 1 ) 


where ^defect is the total energy of the MX2 supercell 
with defects and Fpristine is the total energy of the MX2 
supercell without any defect, ni denotes the number of 
i element (M or X atoms) added to or removed (with a 
negative sign) from the pristine system to create defects. 
jj^i is the chemical potential of the element i. 

To compute the chemical potential, we have used the 
formula {Im + 2 x jix = EMX2 ? where {Im and jix are 
the chemical potentials for M and X respectively. The 
total energy per formula unit of the pristine monolayer 
MX2 is denoted as /xmAs- For M rich environment, 
k'M,max = where is the total energy of a bee 

M metal per atom. Thus, iix.min can be computed as 
Fa, min = (fma2 - For the X rich environ¬ 
ment, Ex,max = where E^^ is the energy per atom 

of the bulk crystal of X. Therefore, Em, min is computed 
as Em, min — EMX2 2 X Ex,max' 


C. Optical properties 


The optical properties are also calculated using the 
VASP code.l^. The optical properties of these materials 
can be described using the complex dielectric constant: 
e{uj) = -h 82(00)• The dielectric function has been 

computed in the momentum representation by obtaining 
the matrix elements between occupied and unoccupied 
eigenstates. The imaginary part of the dielectric function 
can be derived from the following formula® 


82 (^) = lim ^ 2w\^5(ec\^ - e^k - ^) 

q^oo 


c,v,k 


(Fck+eQ,q|Fvk) (Fck+e^gq |Fvk) 5 ( 2 ) 


where the indices c and v denote, respectively, the con¬ 
duction and valence band states. The cell periodic part 
of the wavefunctions at the given k-point is denoted by 
/ick- The real part of the dielectric function is obtained 
by using Kramers-Kronig transformation: 


Using the above parameters we have computed the in¬ 
plane lattice constants and band gaps of different pris¬ 
tine MX2 systems. These computed values along with 
the experimental values have been tabulated in table [I| 
Our calculated lattice constants and band gaps are in 
good agreement with the previous theoretical calcula- 
tionPimsi 


ei(w) = IH —P 

TT 


CXD 



0 


£ 2 (^ 0 ^' . , 

~l2 - 

uo'^ — UO^ d- ITj 


( 3 ) 


where P and 77 denote, respectively, the principal value 
and the complex shift. 
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Figure 1. (Color online) Representative figure of optimized geometries for various defects in MX 2 system as designed on 5x5x1 
supercelL Fig. [^a)-^f) represent the following defects: (a) X-interstitial, (b) X-vacancy, (c) M-vacancy, (d) MX-vacancy, 
(e) XX-vacancy and (f) M-interstitial. The cyan (large) balls denote the X-atoms from the top layer, pink (medium) balls 
denote M-atoms and dark blue (small) balls denote X-atoms from the bottom layer. In Fig. [^a), the purple ball refers the X 
interstitial atom. In Fig.[^f), the red ball indicates to the M interstitial atom. Red arrows in Fig.[^b) indicates the movement 
of M atoms during the relaxation. AL and AH in Fig. [^c) represents the change of in-plane and out-of-plane displacements 
of the X atoms around the M vacancy due to strutural relaxation. Please see section |ra3 for detailed analysis of structural 
relaxations. 


III. RESULTS AND DISCUSSIONS 
A. Geometry &; Energetics 

We have compared six different types of defects for 
various MX2 (M = Mo,W; X = S,Se,Te) systems. The 
defects are as follows - X-vacancy, X-interstitial, M- 
vacancy, M-interstitial, XX-vacancy and MX-vacancy in 
each monolayer MX2 system. To find out the ground 
state structures, we have considered many possible start¬ 
ing geometries to find out the lowest energy structure for 


each of the system. In the following sections, we will first 
discuss the optimized geometries and energetics of point 
defects followed by the double defects. In Fig. we have 
shown the typical lowest energy optimized structures of 
all the different defects that we have considered for the 
present study. 
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Figure 2. (Color online) Formation energies for different types 
of defects in various MX 2 systems under both X-rich and M- 
rich conditions. See subsection IIII Al for more details about 
different systems. 


Point defects 

X - interstitial 

Fig. [^a) represents the energetically most stable X in¬ 
terstitial defect (X^) structure. Our calculation of forma¬ 
tion energy (See Fig. reveals that under both X-rich 
and M-rich conditions, the most stable structure is found 
to have X^ adatom attached to the top of a host X atom 
(X/i). We have considered two other interstitial positions 
for X^ adatom - (i) hexagonal position at the M layer and 
(ii) bridge position between two host X atoms at the top 
layer. For X^ adatom on top of X/^ atom structures, the 
formation energies are ^ 1.1 eV for all MX2 defects for 
X-rich environment. However, for M-rich environment, 
the formation energy drops from ^ 2.4 eV (X = S) to 
^ 1.3 eV (for X = Te). The formation energies for the 
hexagonal position are ^ 6 eV for all MX2 systems, which 
are quite high compared to the X^ adatom attached to 
the top of X/j, atom. Our calculation also shows that the 
above mentioned bridge position is a metastable position 
and after relaxation, the X^ adatom moves to the top of 
Xh adatom. The X^-X/^ bond lengths are 1.94 A, 2.26 A 
and 2.65 A, respectively for X = S, X = Se and X = Te. 


X - vacancy 

The most stable structure of X vacancy is shown in 
Fig. gb). The vacancy is created by removing one X 
atom in the single layer MX2 supercell. The absence 
of one X atom causes the three M atoms to relax to¬ 
wards the vacancy site. M-X bond lengths change by 
0.04 A and 0.02 A respectively for M0X2 and WX2 sys¬ 
tems. Under X-rich environment, the formation energy 
of X vacancy in MoTe2 is 0.21 eV smaller than for MoSe2 
and M0S2. The same for WTe2 is 0.44 eV and 0.33 eV 
smaller than the WSe2 and WS2 respectively. However, 
formation energy of X vacancy in MTe2 is higher than 
MSe2 and MS2 under M-rich environment. 


M - vacancy 

Fig. 0 c) represents the energetically most stable struc¬ 
ture of M vacancy. As the M atom was connected to six 
X atoms, there are six dangling bonds of X atom present 
in the structure affecting the relaxation of these X atoms. 
Table [TT| refers to the length variation due to creation of 
the M vacancy. The analysis of length variation shows 
that while the S and Se atoms relaxed outwards from the 
vacancy center, the Te atoms relaxed inwards to the va¬ 
cancy center. Also, the vertical height between the Te 
atoms from the top and bottom layer reduces by a large 
amount (^ 0.5 A) compared to the S and Se atoms. The 
reason behind these geometry changes is that the S and 
Se are more electronegative compared to Te. Hence the 
dangling S/Se atoms repel each other strongly leading to 
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Table II. Bond length variation due to the creation of M 
vacancy with reference to pristine system. AL is the in-plane 
X-X length variation at the vacancy site and AH is the 
vertical X-X length variation between the two X atoms from 
bottom and top layers at the vacancy site. Am-x is the 
change in the bond length between the X atom with dangling 
bonds and the M atom. 


Length (A) 

M 0 S 2 MoSe 2 MoTe 2 

WS 2 WSe 2 WTe 2 

AL 

0 

0 

0 

0 

0 

1 

0 

to 

0.11 0.05 -0.09 

AH 

-0.02 -0.09 -0.56 

-0.02 -0.08 -0.53 

Am-x 

-0.05 -0.02 0.01 

-0.05 -0.03 0.01 


the outward relaxation from the vacancy center. For Te 
atoms, although they repel each other, there is a bigger 
void to fill in due to the M vacancy. Hence for MTe2 sys¬ 
tem, we observe an inward relaxation. These geometry 
changes directly influence the vacancy formation ener¬ 
gies, where MTe2 has the lowest formation energy com¬ 
pared to MS2 and MSe2 under both X-rich and M-rich 
environment. 


M - interstitial 



Figure 3. (Color online) Close up view for W interstitial de¬ 
fects in WTe 2 . (a) slanted top view, (b) side view. The cyan 
(large) balls denote the X-atoms from the top layer, pink 
(medium) balls denote M-atoms and dark blue (small) balls 
denote X-atoms from the bottom layer. The red ball shows 
the W interstitial atom. 


Table III. M-M^ bond length (in A)variation in M interstitial 
defect structure in MX 2 systems. 


M 0 S 2 MoSe 2 MoTe 2 

WS 2 WSe 2 

2.11 2.06 1.97 

2.24 2.23 


Finally, we show the ground state geometry for M in¬ 
terstitial in Fig. [^f). We have considered two different 


interstitial positions of M atom. In one case, we have in¬ 
serted the M atom at the bridge position in between two 
X atoms and in the other case, the M atom was inserted 
in the split interstitial position along the c direction. Our 
calculated formation energy shows that the split intersti¬ 
tial position is energetically most favorable except WTe2. 
Table m shows the bond length variations for M-M^ for 
the split interstitial position of different MX2 systems. 
Fig-i shows the top and side views of energetically fa¬ 
vorable structure of W interstitial defect in WTe2, which 
is quite different from the structure described above. In 
this case, the interstitial atom resides at the hexagonal 
position in the M layer forming a distorted hexagon. 

Double defects 

MX - vacancy 


Table IV. Bond length variation due to the creation of a MX 
vacancy with reference to pristine system. ALi = change of 
bond length between X2 and X3 atoms. AL 2 = change of 
bond length between X4 and X5. Am-xi = change in the 
bond length between XI and M atoms. 


Length (A) 

M 0 S 2 MoSe 2 MoTe 2 

WS 2 WSe 2 WTe 2 

ALi 

0.22 0.08 -0.35 

0.23 0.16 0.08 

AL 2 

-0.05 -0.26 -0.49 

-0.02 -0.15 0.08 

Am-xi 

-0.12 -0.12 -0.13 

-0.12 -0.12 -0.13 



Figure 4. (Color online) Close up view of optimized geometry 
for WTe vacancy in WTe 2 . (a) slanted top view, (b) side view. 
The cyan (large) balls denote the X-atoms from the top layer, 
pink (medium) balls denote M-atoms and dark blue (small) 
balls denote X-atoms from the bottom layer. 


For a MX di-vacancy (created by removing adjacent 
M and X atoms), the most stable structure is shown in 
Fig. [^d). Due to the removal of M and X atoms, five 
neighboring X atoms around the vacancy site relax dur¬ 
ing the optimization process. Table [IV| shows the amount 
of relaxation of these X atoms with dangling bonds. Af¬ 
ter geometry optimization, XI atom in general moves 
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towards the M layer reducing the M-Xl bond length for 
all the systems. For MX2 (X = S, Se), the X2 and X 3 
atoms relax outward from the vacancy site while X 4 and 
X 5 atoms relax inward mainly because of the missing X 
atom in that layer. The reason for these relaxation be¬ 
havior is similar to the case of M vacancy. For MoTe2, 
the dangling X atom relaxes inward. However, for WTe2 
system they relax slightly outward because of the fact 
that XI atom moves to M layer and forms an isosceles 
triangle with two M atoms with bond length 2.6 A. Fig.|^ 
shows the slanted top and side views of this particular ge¬ 
ometry. The formation energy of MX vacancy is lower 
for MTe2 for both X and M rich environment. 


XX - vacancy 

The most stable structure of XX di-vacancy is shown 
in Fig. [^e). We have considered two different configura¬ 
tions for XX vacancy. In one configuration, we have re¬ 
moved two X atoms (with same x and y coordinates) from 
the top and bottom layers (vac-XXl) and in the other 
configuration, two nearest X atoms have been removed 
from the same layer (vac-XX 2 ). The calculation of for¬ 
mation energies indicates that the vac-XXl structure is 
energetically more favorable compared to vac-XX 2 struc¬ 
ture. Here we find that the M atoms relax inward to the 
vacancy site and form an equilateral triangle with bond 
length ^ 2.8 A. 

Table V. Relative defect formation energies (in eV) for XX- 
and MX-vacancies for different distances between vacancies 
(upper panel of the table). The corresponding distances 
(in A) between the vacancies (before relaxation) are also 
mentioned in the bottom panel of the table, where M stands 
for Mo and W. 



XXI1 

XXI2 

XXI3 

XX2i 

XX22 

MXi 

MX2 

M0S2 

0.0 

0.13 

0.04 

0.00 

0.05 

0.0 

2.02 

MoSe2 

0.0 

0.45 

0.32 

0.24 

0.35 

0.0 

1.88 

MoTe2 

0.0 

0.99 

0.83 

0.64 

0.86 

0.0 

1.59 

WS2 

0.0 

0.31 

0.18 

0.2 

0.21 

0.0 

1.89 

WSe2 

0.0 

0.68 

0.50 

0.47 

0.54 

0.0 

1.95 

WTe2 

0.0 

1.32 

1.03 

0.96 

1.08 

0.0 

1.84 


XXI1 

XXI2 

XXI3 

XX2i 

XX22 

MXi MX2 

MS2 

3.12 

4.46 

6.33 

3.18 

5.59 

2.41 

3.99 

MSe2 

3.34 

4.71 

6.65 

3.32 

5.75 

2.54 

4.18 

MTe2 

3.60 

5.07 

7.14 

3.56 

6.17 

2.73 

4.49 


To investigate the dependence of distance on the for¬ 
mation energies of XX- and MX-vacancies, we have cal¬ 
culated defect formation energies for various distances 
between a pair of defects. In table |V| relative energies as 


a function of distances are tabulated. XXI and XX 2 va¬ 
cancies denote two X atoms removed from different layers 
and same layers of X atom, respectively, from a pristine 
MX2 system. XXli, XXI2 and XXI3 are three different 
XXI configurations considered with three different dis¬ 
tances between X atoms where XXli corresponds to the 
nearest neighbor configuration. In a same manner, XX 2 i 
and XX22 are two different XX 2 configurations consid¬ 
ered. For MX vacancy also we have considered two differ¬ 
ent distances - MXi and MX2 where MXi is the nearest 
neighbor. Our results clearly show that the double va¬ 
cancies always prefer to form when they are within the 
nearest neighbor distance signifying the natural occur¬ 
rence of correlated vacancies. A clear trend of increasing 
formation energies while one goes down in the periodic 
table from S to Te is observed for both Mo and W in case 
of anion vacancies. The trend is not absolutely clear for 
M-X vacancies. 


B. Equilibrium defect concentration 

We have also calculated equilibrium defect concentra¬ 
tions for all the defects. The equilibrium defect concen¬ 
tration, Ceg, can be calculated using the occupation prob¬ 
ability as shown below, 

Ce, = Ne-^CksT (4) 

where N is the concentration of possible defect sites, Ef 
is the formation energy of the defect and T is the tem¬ 
perature (lOOOK - 1200 K) during the crystal growth. 
In Tables [Vl] and 1^ we have tabulated the range of 
equilibrium defect concentrations for a range of growth 
temperatures (lOOOK - 1200 K) for X-rich and M-rich 
conditions respectively. This range of temperature is cho¬ 
sen as most of the MX2 monolayer structures are synthe¬ 
sized experimentally in this range.^^^M^. Our calculation 
clearly shows that the presence of Xint defect under the 
X-rich condition and X^ac defect under the M-rich con¬ 
dition during the crystal growth is very much likely. How¬ 
ever the probability of forming Mint defect under X-rich 
condition and the formation of Myac defect under M-rich 
condition is very low. 


C. Electronic properties 

In this section, we will discuss the electronic properties 
of the defected MX2 structures by analyzing total density 
of states (DOS), site and mi projected densities of states 
(PDOS) and partial charge densities. 

In Fig. we have shown an energy level diagram for all 
the different defect states in all MX2 systems constructed 
from our calculated data. In each subplot, the columns 
(a)-(f) represent the (a) X-interstitial, (b) X-vacancy, 
(c) M-vacancy, (d) MX-vacancy, (e) XX-vacancy and 
(f) M-interstitial defect respectively. The orange long 
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Table VI. Equilibrium defect concentration, Ceg, in cm ^ for X-rich condition for growth temperature in the range of 
1000K-1200K 




X, 






MX^ 


XX, 


M 0 S 2 

MoSe2 

MoTe2 

WS 2 

WSe 2 

WTe2 


9.2x10^-7.2x10 

1.0x10^-1.1x10 

5.1x10^-4.5x10 

1.2x10^-1.4x10 

2.2x10^-3.3x10 

3.0x10^-2.9x10 


10 

10 

10 

10 

9 

10 


96.7- 1.6x10"^ 5.1x10“^®-!.8x10“^® 

1.1 X 10^-1.8 X 10"^ 5.4x 10“^’'-8.9x 10“^^ 
1.1x10^-1.2x10® 7.6x10““- 1.2x10“® 

13.7- 3.2x10® 1.7xl0“^®-7.5xl0“®® 

48.6-9.2x10® 1.3x10“^^-!.8x10“®® 

2.2xl0®-2.2xl0® 0.1-47.6 


8.4xl0“®-4.1xl0“"‘ 

5.3xl0“"‘-0.6 

67.1- 1.1x10® 
8.9x10“®-!.2x10“^ 
1.5xl0“^-9.78 

32.2- 5.8x10® 


6.3x10“®®-!.1x10“® 

1.6xl0“®'-7.9xl0“® 

9.9x10“®-!.14 

1.2xl0“®°-2.0xl0“® 

1.2xl0“®’-6.4xl0“® 

0.3-1.1x10^ 


2.7xl0“®^-7.3xl0“® 
7.5x10“®®-!.2x10“® 
1.8xl0“®-5.2xl0“® 
3.1x10“®®-!.2x10“® 
1.4x10“®°-!.9x10“® 
1.3xl0“®-0.2 


Table VII. Equilibrium defect concentration, Ceq, in cm ^ for M-rich condition for growth temperature in the range of 
1000K-1200K 



'^int 

^vac 

Mint 

M,a,c 

MX,ac 

XX,ac 

M 0 S 2 

MoSe2 

MoTe2 

WS 2 

WSe2 

WTe2 

1.8xl0®-1.9xl0® 

5.8xl0®-5.0xl0® 

4.8xl0®'-9.1xl0® 

6.0x10^-7.5x10® 

I.lxl0®-8.4xl0® 

1.0xl0°-1.2xl0®° 

4.8xl0®-6.2xl0° 

1.8xl0®'-4.1xl0® 

1.2xl0®-6.0xl0® 

2.8xl0®'-5.8xl0® 

9.8xl0®-3.6xl0®' 

6.4xl0®-5.4xl0® 

1.2xl0“®-2.6xl0“^ 

1.6xl0“®-4.7xl0“® 

8.8xl0“®'-2.9xl0“® 

7.3xl0“®®-2.5xl0“® 

5.2xl0“®®-2.7xl0“° 

0.9-2.9x10^ 

3.4x10“2®-2.9x10“®® 

1.8x10“®®-!.1x10“° 

5.8xl0“®-4.4 

2.1xl0“®®-6.0xl0“®® 

3.7xl0“®®-6.5xl0“®' 

3.7-9.6x10^ 

1.3xl0“®®'-3.0xl0“®^ 

9.2xl0“®®-3.4xl0“® 

9.2x10“®-2.3x10“2 

5.9x10“®®'-!.1x10“®® 

6.1x10“®M.6x10“®' 
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Figure 5. (Color online) Energy level diagram for different defects in all the MX 2 systems. In each subplot, colums (a)-(f) refer 
to X-interstitial, X-vacancy, M-vacancy, MX-vacancy, XX-vacancy and M-interstitial respectively. The orange (long dashed 
in print) line shows the position of VBM and CBM for pristine MX 2 . Green solid and dashed lines denote the occupied and 
unoccupied defect states respectively. 
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Figure 6. (Color online) Total and partial densities of states 
for X interstitial in M 0 S 2 and MoTe 2 . Fig.j^a) shows partial 
charge density of MoTe 2 showing electronic orbitals for VBM. 
Fig.^b) shows the same for CBM. Black (solid) line denotes 
the total density of states. Green (checker shaded) region 
represents the defect states due to Mo-d orbitals. Orange 
(solid shaded) region shows the defect states due to S,Te-p 
states. 


dashed line denote the position of valence band maxima 
(VBM) and conduction band minima (CBM) of the pris¬ 
tine MX 2 systems. The short green solid and dashed 
lines denote the position of defect states for both occu¬ 
pied and unoccupied states respectively. From the analy¬ 
sis of energy level diagram, it can be seen that the defect 
states mainly appear in the gap region of the pristine sys¬ 
tem. Further analysis of DOS for all the defects in all the 
MX 2 systems reveals that for a specific defect, the quali¬ 
tative behavior is very much similar between M 0 S 2 /WS 2 , 
MoSe 2 /WSe 2 and MoTe 2 /WTe 2 . Furthermore, our anal¬ 
ysis shows that the qualitative electronic properties due 
the presence of defects in MoS 2 /MoSe 2 are quite similar. 
Hence we have shown the DOS, PD OS analysis for only 
M 0 S 2 MoTe 2 systems. 

The analysis of density of states for X interstitial de¬ 
fects for all the systems shows the same qualitative be¬ 
havior. In Fig. I^c), we have shown the total densities of 
states for two systems - M 0 S 2 and MoTe 2 . The general 
observation is that the defect states are found near the 
band edges. These states merge with both valence and 
conduction bands. To find the orbital character of these 
defect states, we have carried out the analysis of site pro¬ 
jected density of states (PDOS) and partial charge den- 
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Figure 7. (Color online) Total and partial DOS for X va¬ 
cancy in M 0 S 2 and MoTe 2 . Fig. [3a) -(c) show partial charge 
densities of different defect peaks in the band gap region as 
indicated in the DOS. See Fig. |^for other notations. 
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Figure 8. (Color online) Total and projected DOS for M va¬ 
cancy in M 0 S 2 and MoTe 2 . Fig.[^a)-(d) show partial charge 
densities for defect peaks in the band gap region as noted in 
the DOS. See Fig. [^for other notations. 
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Figure 9. (Color online) Total and partial densities of states 
for MX vacancy in M 0 S 2 and MoTe 2 . Fig. [^a)-(g) show 
partial charge densities for defect peaks in the band gap region 
as noted in the DOS. See Fig. |^for other notations. 


sity analysis. Figs.j^a) and[^b) show the partial charge 
density for valence band maximum (VBM) and conduc¬ 
tion band minimum (CBM) respectively. Our analysis 
shows that the states near the VBM mainly arise from 
atom. These states are the doubly degenerately oc¬ 
cupied Px and Py orbitals. There are also very small 
contributions of in-plane M - dxy and dx 2 -y 2 states ap¬ 
pearing in the VBM. The states near the CBM originate 
mainly from the M - orbitals with some contribution 
from the unoccupied Pz orbital of X^ atom. 

Fig. shows density of states for X vacancy defects for 
M 0 S 2 and MoTe 2 . For X vacancy defects, the features of 
density of states in the valence band side are identical for 
all the MX 2 systems. In these systems, the defect state 
appears at the valence band edge and inside the band gap 
towards the conduction band side. The defect states can 
be further classified as singlet state (near valence band) 
and doubly degenerate state (the empty localized state in 
the gap) due to the trigonal symmetry of the X vacancy. 
The analysis of on-site projected DOS reveals that these 
states are formed from the combination of diagonal out of 
plane d orbitals {dxz and dyz) and in-plane d orbitals {dxy 
and dx 2 -y 2 ) of the M atoms near the defected sites and 
very small amount of in-plane p orbitals (mainly Px) of 




- 2-1012 
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Figure 10. (Color online) Comparison of total density states 
and defect states for XX vacancy in M 0 S 2 and MoTe 2 . 
Fig. p^a)-(c) show partial charge density for defect peaks 
in the band gap region as noted in the DOS. See Fig. [^for 
other notations. 
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Figure 11. (Color online) Comparison of total density states 
and defect states for M interstitial in M 0 S 2 and MoTe 2 . 
Fig. ll^a)-(c) show partial charge density for defect peaks 


in the band gap region as noted in the DOS. See Fig. [^for 


other notations. 
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second nearest neighbor X atoms. The orbital characters 
of VBM and CBM are respectively contributed by the 
dx 2 _y 2 and (i ^2 orbitals of the M atom. 

The analysis of PDOS and partial charge density shows 
that the defect states due to the M vacancy appear inside 
the gap region. The density of states for M vacancy in 
M 0 S 2 and MoTe 2 systems are shown in Fig. In this 
case, singlet and doublet defect states can be seen due 
to the trigonal and mirror symmetry of M vacancy and 
two layers of X atoms, respectively. The VBM contains 
both the Pz orbital (significant contribution) of the six X 
atoms near the vacancy site and (1^2 _y 2 orbitals of the M 
atoms (see Fig.j^a)). The defect peak (marked as (b) in 
DOS) is a doublet state and originates mainly from the 
dx 2 _y 2 and dz 2 orbitals of the six M atoms (second near 
neighbor from the vacancy site). The other defect peak 
(marked as (c) in the DOS) originates from the Pz orbital 
of the six X atoms near the vacancy site along with very 
small dxz orbital contribution from the M atoms. The 
CBM mainly contains the dx 2 -y 2 and dz^ orbitals of the 
M atoms. All these states are formed due to the local 
perturbation created in the geometry due to the absence 
of one M atom. 

In Fig. we have shown the MX-vacancy density of 
states comparison for M 0 S 2 MoTe 2 as well as the 
partial charge density for the defect states, VBM and 
CBM. For MX vacancy, combined defect states due to 
M-vacancy and X-vacancy can be observed. As seen for 
the figure (see Fig. [^, the defect states appear mainly 
in the gap region. This defect has a strong influence 
on conduction band of MX 2 . The states marked as (a) 
and (b) in the DOS) are mainly originating from the 
dxy and dx 2 -y 2 orbital of the M atoms near the vacancy 
site. A mix orbital character of dxy^ dz^ form the M 
atom and in-plane p orbital from the dangling X atom 
can be observed in defect peak (c). However, the major 
contribution to the defect peak (d) is coming from the Pz 
orbitals of the dangling X atom with a minor contribution 
from dyz orbital of M atom. The defect peak (e) is again 
a mixture of dx 2 -y 2 , dz 2 of M atom and Pz of the X atom. 
However, the defect peak (f) is mainly of d character with 
dyz and dx 2 -y 2 are being the dominant orbitals. The 
CBM mainly consists of dz 2 orbital of M atoms. 

For the XX double vacancy also the defect states can 
be seen in the gap region. Fig. pT| shows the comparison 
of densities of states for M 0 S 2 and MoTe 2 and the partial 
charge densities of MoTe 2 system for this defect. The 
analysis of partial charge density and PDOS shows that 
these defect states are the combined effect of two single X 
vacancies. Consequently, the mi characteristic features of 
the defect peaks are very much similar to the X-vacancy 
defects as discussed previously. 

The density of states for Mo interstitial in M 0 S 2 
MoTe 2 and the partial charge density for the MoTe 2 have 
been shown in Fig. The general trend in the density 
of states for Mo interstitial defect in MX 2 system is again 
very similar for all the systems. The defect levels are gen¬ 
erated inside the band gap. The occupied states are one 


doublet and one singlet state and the unoccupied state 
is the doublet state. The contribution to the occupied 
defect states are mainly coming from the dz 2 orbital of 
the interstitial M atom and the M atom attached with it. 
There is also a smaller contribution from the neighboring 
M atom’s dx 2 -y 2 orbital. However, the contribution to 
the unoccupied defect state is manly from the defect site 
M atom’s dxy and dx 2 -y 2 orbitals. 


D. Optical properties 



Energy (eV) Energy (eV) Energy (eV) 


Figure 12. (Color online) Optical properties of pure MX 2 sys¬ 
tem. Blue and orange lines denote In-plane and out-of-plane 
components of the imaginary parts of the dielectric functions 
(£ 2 ). 


In addition to the geometric and electronic structures, 
we have also studied the influence of defects in the optical 
properties of defected MX 2 by means of frequency depen¬ 
dent dielectric functions. In the following paragraphs we 
will briefly discuss the comparison of optical properties 
for the pristine system followed by the modification of 
optical properties due to the presence of defects. 

The comparison of optical properties for all the pris¬ 
tine MX 2 systems have been shown in Fig. We have 
plotted the imaginary parts of the dielectric functions 
with their in-plane {xx and yy) and out of plane (zz) 
components as a function of energy. The xx and yy com¬ 
ponents have the same absorption spectra for all the sys¬ 
tems and their absorption is stronger than the 2 : 2 ; com¬ 
ponent. From the figure, it can be seen that the general 
features of the optical spectra for same X elements are 
quite similar in nature with a step like peak followed by a 
dominant peak, valley and subsequent peaks. This step 
like peak is the characteristic of a two-dimensional sys¬ 
tem. From the analysis of the spectra, it can be seen that 
the prominent optical peak shifts towards lower energy 
and breaks into more than one peak as we move towards 
the heavier X element (see fig. 12). The step like peak 
also moves towards lower energy with heavier X element, 
signifying smaller band gaps. 
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Figure 13. (Color online) Comparison of imaginary parts of 
the dielectric functions (£ 2 ) between M 0 S 2 , MoSe 2 and MoTe 2 
for M-vacancy, X-vacancy and MX-vacancy. Blue and or¬ 
ange lines denote in-plane contribution of £2 for defected and 
pure systems respectively. Purple and green lines denote out- 
of-plane contribution of £2 for defected and pure systems re¬ 
spectively. The shaded region represents the prominent con¬ 
tributions of defect related peaks compare to the pristine sys¬ 
tem. 


M-Vacancy X-Vacancy MX-Vacancy 



Energy (eV) Energy (eV) Energy (eV) 


Figure 14. (Color online) Comparison of optical properties 
between WS 2 , WSe 2 and WTe 2 for M-vacancy, X-vacancy 
and MX-vacancy. See Fig. for other notations. 


Introduction of defects in the pristine MX 2 system can 
affect the optical properties also. From our calculations, 
we have found out that most notable change in the optical 
properties due to the presence of defect occurs for three 
following defects - i) M-vacancy, ii) X-vacancy and iii) 
MX-vacancy. Therefore, in the following section we have 
discussed aforementioned defects only. 

In Fig. Poland Fig.we have shown the comparison 
of optical properties for M-, X- and MX-vacancies for 
M 0 X 2 and WX 2 respectively. As mentioned before, here 
also the xx and yy components have the same absorp¬ 
tion spectra for all the systems and their absorption is 
stronger than the 2 : 2 : component. 


For M-vacancy in M 0 S 2 , defect peaks appear just 
below and above the Fermi level (see fig. |^, which have 
mainly p and d characters and are responsible for the 
electronic transition. This transition occurs at ^ 1.0 eV. 
In MoSe 2 , the peak is broader as there are two possi¬ 
ble transitions possible with a very small energy inter¬ 
val. However for M-vacancy in MoT^, the defect peak 
appears at the Fermi energy (see fig. giving rise to a 
sharp optical transition at a very low energy (see fig. 13). 


For X-vacancy, the defect state appears very close to 
the CBM of pristine system (see fig. [7). Thus the optical 
spectra with X-vacancy in different M 0 X 2 systems do not 
change significantly from the pristine one. In the case of 
MX-vacancy also a number of defect peaks appear near 
the Fermi energy for all M 0 X 2 systems (See fig. which 
give rise to the electronic transition at ^ 0.5 eV with a 
relatively broader peak distribution than M-vacancy. 


The analysis of DOS for WX 2 systems (figure not 
shown) show similar characteristics as of M 0 X 2 . This 
leads to a very much similar optical properties as shown 
in Fig. as compared to M 0 X 2 (Fig, [l^. 

The calculated dielectric functions elicit that the ab¬ 
sorption spectrum appears in the visible region (3 eV-1.7 
eV (400-700 nm)) for MX 2 systems. It is well known that 
the material, which has absorption in the visible region 
is suitable for photocatalysis using sunlight. Therefore, 
MX 2 having suitable defects can be used for the photo¬ 
catalysis. Tongay et al ^ in their study have investi¬ 
gated effects of anion vacancy on the photoluminescence 
of M 0 S 2 , MoSe 2 and WS 2 where this vacancy is identified 
by a new peak appearing below the photoluminescence 
peak. Moreover, the intensity of the photoluminescence 
peak enhances with increasing defect density. As we dis¬ 
cussed above, the absorption spectra of both M-vacancy 
and MX-vacancy defected MX 2 have new peaks at ^ 1.0 
eV and ^ 0.5 eV, respectively, which arise due to the 
defect states. This feature is very much similar to that 
observed with anion vacancy. Nowadays, it is possible 
to create point defects in 2D materials (as demonstrated 
in graphene) using ion-irradiation technique. One may 
speculate that the creation of suitable defects in MX 2 
may give rise to desired optical transitions suitable for 
light emitting diodes. 
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IV. CONCLUSIONS 


We have performed a systematic study of native de¬ 
fects, viz. X vacancy, X interstitial, M vacancy, M in¬ 
terstitial, MX vacancy and XX vacancy on a single layer 
of MX 2 system (M = Mo, W; X = S, Se, Te) using first 
principles density functional theory. It has been found 
that under X-rich condition, X interstitial defect has the 
lowest formation energy for all the systems under inves¬ 
tigation with equilibrium defect concentration ^ xlO^ 
cm“^ at the growth temperature range of 1000-1200K. 
For M-rich environment, X vacancy has the lowest for¬ 
mation energy except MTe 2 systems, where X interstitial 
defect has the lowest formation energy. Metal atom and 
double vacancy related defects are quite high in forma¬ 
tion energy thus showing almost zero equilibrium defect 
concentration at the studied growth temperatures. Our 
calculations reveal the position of defect states in the 
band gap along with the orbitals contributing to those 


states. Finally, our calculated optical properties indi¬ 
cate prominent signatures in the in-plane component of 
the imaginary part of the dielectric functions. One may 
speculate that suitably designed defected MX 2 systems 
can be promising materials for light emitting devices. 

ACKNOWLEDGEMENT 

OE acknowledges KAW foundation for financial sup¬ 
port. In addition, BS acknowledges Carl Tryggers Stif- 
telse, Swedish Research Council and KOF initiative of 
Uppsala University for financial support. SH also ac¬ 
knowledges valuable discussions with Carmine Autieri. 
We are grateful to NSC under Swedish National Infras¬ 
tructure for Computing (SNIC) and the PRACE project 
resource Cy-Tera supercomputer based in Cyprus at the 
Computation-based Science and Technology Research 
Center (CaSToRC) and Zeus supercomputer based in 
Poland at the academic computer center Cyfronet. Struc¬ 
tural figures were generated using VMEE^. 


* Soumyaj yoti. Haidar@physics .uu.se 

^ Current address: Nepal Academy of Science and Technol¬ 
ogy, Khumaltar, Lalitpur, Nepal. 

^ Corresponding author: Biplab.Sanyal@physics.uu.se 
^ K. S. Novoselov, A. K. Geim, S. V. Morozov, 

D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grig¬ 

orieva, and A. A. Firsov, Science 306, 666 (2004), 
http://www.sciencemag.org/content /306/5696/666.full.pdf 
^ A. O’Hare, F. V. Kusmartsev, and K. I. Kugel, 

Nano Letters 12 , 1045 (2012), pMID: 22236130, 

http://dx.doi.org/I0.I02I/nl204283q 
^ H. Okamoto, Y. Kumai, Y. Sugiyama, T. Mitsuoka, 

K. Nakanishi, T. Ohta, H. Nozaki, S. Yamaguchi, 
S. Shirai, and H. Nakano, Journal of the American 
Chemical Society 132, 2710 (2010), pMID: 20I2I277, 
http://dx.doi.org/I0.I02I/ja908827z 
^ Y. Sugiyama, H. Okamoto, T. Mitsuoka, T. Morikawa, 
K. Nakanishi, T. Ohta, and H. Nakano, Journal 
of the American Chemical Society 132, 5946 (2010), 
http://dx.doi.Org/I0.I02I/jaI009I9d 
^ X. Yang and J. Ni, Phys. Rev. B 72, 195426 (2005) 

® S. Rachel and M. Ezawa, Phys. Rev. B 89, 195303 (2014), 

^ L. Li, S.-z. Lu, J. Pan, Z. Qin, Y.-q. Wang, Y. Wang, G.-y. 
Cao, S. Du, and H.-J. Gao, Advanced Materials 26, 4820 
(2014) 

^ Y. Cai, C.-P. Chuu, C. M. Wei, and M. Y. Chou, Phys. 
Rev. B 88, 245408 (2013), 

Y. Xu, B. Yan, H.-J. Zhang, J. Wang, G. Xu, P. Tang, 
W. Duan, and S.-C. Zhang, Phys. Rev. Lett. Ill, 136804 
(2013) 

B. van den Broek, M. Houssa, E. Scalise, G. Pourtois, 
V. V. Afanas‘ev, and A. Stesmans, 2D Materials 1 , 021004 
(2014) 

H. Liu, A. T. Neal, Z. Zhu, Z. Luo, X. Xu, D. Tomanek, 
and P. D. Ye, ACS Nano 8, 4033 (2014), pMID: 24655084. 


S. Das, W. Zhang, M. Demarteau, A. Hoffmann, M. Dubey, 
and A. Roelofs, Nano Letters 14, 5733 (2014) 

M. M. Ugeda, A. J. Bradley, S.-F. Shi, F. H. da Jornada, 

Y. Zhang, D. Y. Qiu, W. Ruan, S.-K. Mo, Z. Hussain, 

Z. -X. Shen, F. Wang, S. G. Louie, and M. F. Crommie, 
Nature Materials 13, 1091 (2014). 

K. F. Mak, C. Lee, J. Hone, J. Shan, and T. F. Heinz, 
Phys. Rev. Lett. 105, 136805 (2010) 

H. Zeng, J. Dai, W. Yao, D. Xiao, and X. Cui, Nature 
Nanotechnology 7, 490 (2012). 

W. Zhou, X. Zou, S. Najmaei, Z. Liu, Y. Shi, J. Kong, 
J. Lou, P. M. Ajayan, B. 1. Yakobson, and J.-C. 
Idrobo, Nano Letters 13, 2615 (2013), pMID: 23659662, 
http://dx.doi.Org/10.1021/nl4007479 

J. -Y. Noh, H. Kim, and Y.-S. Kim, Phys. Rev. B 89, 
205417 (2014) 

Y.-H. Lee, X.-Q. Zhang, W. Zhang, M.-T. Chang, C.-T. 
Lin, K.-D. Chang, Y.-C. Yu, J. T.-W. Wang, C.-S. Chang, 

L. -J. Li, and T.-W. Lin, Advanced Materials 24, 2320 

( 2012 ) 

B. Radisavljevic, A. Radenovic, J. Brivio, V. Giacometti, 
and A. Kis, Nature Nanotechnology 6, 147 (2011). 

K. Q. Dang and D. E. Spearot, Journal of Applied Physics 
116, 013508 (2014) 

S. Tongay, J. Suh, C. Ataca, W. Fan, A. Luce, J. S. Kang, 
J. Liu, C. Ko, R. Raghunathanan, J. Zhou, F. Ogletree, 
J. Li, J. C. Grossman, and J. Wu, Sci. Rep. 3, 2657 (2013) 

L. ping Feng, J. Su, and Z. tang Liu, Journal of Alloys 
and Compounds 613, 122 (2014) 

M. Ghorbani-Asl, A. N. Enyashin, A. Kuc, G. Seifert, and 

T. Heine, Physical Review B 88, 245440 (2013). 

H.-P. Komsa, S. Kurasch, O. Lehtinen, U. Kaiser, and 
A. V. Krasheninnikov, Physical Review B 88, 035301 
(2013). 

H.-P. Komsa and A. V. Krasheninnikov, Physical Review 
B 91, 125304 (2015). 




13 


S. KC, R. C. Longo, R. Addou, R. M. Wallace, and 
K. Cho, , 1 (2014). 

C. Ataca and S. Ciraci, The Journal of Physical Chemistry 
C 115, 13303 (2011). 

Y.-C. Lin, T. Bjorkman, H.-P. Komsa, P.-Y. Teng, C.-H. 
Yeh, F.-S. Huang, K.-H. Lin, J. Jadczak, Y.-S. Huang, P.- 

W. Chiu, A. V. Krasheninnikov, and K. Suenaga, Nat 
Commun 6 (2015). 

G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996) 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996) 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 78, 1396 (1997) 

T. Boker, R. Severin, A. Muller, C. Janowitz, R. Manzke, 

D. VoB, P. Kruger, A. Mazur, and J. Pollmann, Phys. Rev. 
B 64, 235305 (2001) 

W. Schutte, J. D. Boer, and E. Jellinek, Journal of Solid 
State Chemistry 70, 207 (1987) 

K. E. Mak, K. He, C. Lee, G. H. Lee, J. Hone, T. E. Heinz, 
and J. Shan, Nature Materials 12 , 207 (2012) 

R. Soklaski, Y. Liang, and L. Yang, Applied Physics Let¬ 
ters 104, 193110 (2014) 

Y. Zhang, T.-R. Chang, B. Zhou, Y.-T. Cui, H. Yan, 

Z. Liu, E. Schmitt, J. Lee, R. Moore, Y. Chen, H. Lin, 

H. -T. Jeng, S.-K. Mo, Z. Hussain, A. Bansil, and Z.-X. 
Shen, Nature nanotechnology 9, 111 (2014) 

M. M. Ugeda, A. J. Bradley, S.-E. Shi, E. H. da Jornada, 

Y. Zhang, D. Y. Qiu, W. Ruan, S.-K. Mo, Z. Hussain, 

Z. -X. Shen, E. Wang, S. G. Louie, and M. E. Crommie, 
Nature Materials 13, 1091 (2014), 1404.2331 

C. Ruppert, O. B. Aslan, and T. E. Heinz, Nano letters 
14, 6231 (2014). 

Y. Wang, C. Cong, W. Yang, J. Shang, N. Peimyoo, 
Y. Chen, J. Kang, J. Wang, W. Huang, and T. Yu, Nano 
Research 8, 2562 (2015) 

S. B. Desai, G. Seol, J. S. Kang, H. Pang, 
C. Battaglia, R. Kapadia, J. W. Ager, J. Guo, 


and A. Javey, Nano Letters 14, 4592 (2014), 

http://dx.doi.org/10.1021/nl501638a 
Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman, 
and M. S. Strano, Nature Publishing Group 7, 699 (2012). 

A. Ramasubramaniam, Physical Review B 86, 115409 

( 2012 ). 

J. Kang, S. Tongay, J. Zhou, J. Li, and J. Wu, Applied 
Physics Letters 102 , 012111 (2013). 

M. Gajdos, K. Hummer, G. Kresse, J. Eurthmiiller, and 

E. Bechstedt, Phys. Rev. B 73, 045112 (2006) 

X. Ling, Y.-H. Lee, Y. Lin, W. Pang, L. Yu, M. S. Dressel- 
haus, and J. Kong, Nano Letters 14, 464 (2014), pMID: 
24475747, http://dx.doi.org/10.1021/nl4033704 
S. Wang, Y. Rong, Y. Pan, M. Pacios, H. Bhaskaran, 

K. He, and J. H. Warner, Chemistry of Materials 26, 
6371 (2014), http://dx.doi.org/10.1021/cm5025662, 

X. Wang, Y. Gong, G. Shi, W. L. Chow, K. Keyshar, 
G. Ye, R. Vajtai, J. Lou, Z. Liu, E. Ringe, B. K. Tay, and 
P. M. Ajayan, ACS Nano 8, 5125 (2014), pMID: 24680389, 
http://dx.doi.org/10.1021/nn501175k 

N. R. Pradhan, D. Rhodes, S. Peng, Y. Xin, S. Memaran, 

B. -H. Moon, H. Terrones, M. Terrones, and L. Bal- 
icas, ACS Nano 8, 5911 (2014), pMID: 24878323, 
http://dx.doi.Org/10.1021/nn501013c 

Y. Zhang, Y. Zhang, Q. Ji, J. Ju, H. Yuan, J. Shi, T. Gao, 
D. Ma, M. Liu, Y. Chen, X. Song, H. Y. Hwang, Y. Cui, 
and Z. Liu, ACS Nano 7, 8963 (2013), pMID: 24047054, 
http://dx.doi.org/10.1021/nn403454e 

J.-K. Huang, J. Pu, C.-L. Hsu, M.-H. Chiu, Z.-Y. Juang, 

Y.-H. Chang, W.-H. Chang, Y. Iwasa, T. Takenobu, 
and L.-J. Li, ACS Nano 8, 923 (2014), pMID: 24328329, 
http://dx.doi.org/10.1021/nn405719x 
W. Humphrey, A. Dalke, and K. Schulten, Journal of 
Molecular Graphics 14, 33 (1996). 


